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Abstract 

We have considered the properties of anisotropic two-leg ladder models with 
5 = 1/2 or 5 = 1 spins on the rungs, using White's density matrix renor- 
malization group method. We have generalized the method by taking into 
account the symmetries of the model in order to reduce the dimensions of the 
matrix to be diagonalized, thereby making possible to consider more states. 
The boundaries in the parameter space of the extended region, where the 
Haldane phase exists, are estimated. 
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I. INTRODUCTION 



The properties of magnetic systems where the localized moments form a ladder like struc- 
ture |I|-£§, have been intensively studied recently. Most of such materials can be described 
by an isotropic spin Hamiltonian, therefore the theoretical studies were almost exclusively 
done on isotropic models PHH1|| . The question arises as to what effect the anisotropy might 
have. The answer to the question is not trivial. 

It is known that two-leg ladder models with S — 1/2 spins on both legs behave basically 
like an S = 1 spin model if the rung coupling between the spins on the two legs is strong 
enough. It is also known that integer and half-odd integer spin models have essentially 
different phase diagrams, when considered as a function of the anisotropy of the exchange. 
Beside the ferromagnetic, antiferromagnetic and planar phases a new one, the Haldane phase 
|TT| appears for a finite range of anisotropy. Our aim in this paper is to study, how this 
new phase appears in anisotropic ladder models for intermediate values of the interchain 
couplings. 

For this purpose we have determined the low lying part of the energy spectrum of an 
anisotropic Heisenberg ladder model by using the density matrix renormalization group 



(DMRG) method proposed recently by White |I2[. We have generalized the method in such 
a way that together with any spin configuration that is kept after the truncation of the 
Hilbert space, the other configurations related to it by the symmetries of the ladder model 
are also automatically taken into account. This allows to consider more states without 
increasing the size of the matrices and thereby to improve the accuracy. 

The layout of the paper is as follows. In Sec. II we give a short description of the ladder 
model and show what phase diagrams are expected, if the value of the spins on the legs 
is 5 = 1/2 or 5 = 1. The symmetry considerations introduced in the application of the 
DMRG procedure are briefly discussed in Sec III. Sec. IV presents the numerical results. 
Finally Sec V. contains a brief discussion of the results. 
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II. LADDER MODELS AND THEIR EXPECTED PHASE DIAGRAM 



In a two- leg ladder model the spins at rung i will be denoted by <7j and T{. The length 
of the spin can be arbitrary. In most of the calculations we will take 5 = 1/2 spins, but a 
comparison with the S = 1 case will also be considered. 

When the coupling between the legs is neglected, the Hamiltonian of two decoupled 
anisotropic Heisenberg chains is recovered. It can be written in the form, 



^ = E 



-J xy (or+cr i+1 + a i af +1 ) + -].n)a 



z z 
z"i a i+l 



+ E [2 J xy(^ T i+i + r * r m) + J * T i T Ui 



(1) 



The xy part of the exchange is written in terms of the raising and lowering operators, and 
its coupling is allowed to differ from that of the zz part. 

Introducing now the couplings between the chains, different kinds of ladder models can 
be constructed depending of the choice of this coupling. Usually this coupling is assumed 
to act between spins on the same rung only. In an anisotropic model this would mean a 
coupling of the form, 



(2) 



In the limit when this coupling is ferromagnetic and strong, the two spins on the same rung 
couple into a single larger spin. That is the reason why a spin- 1/2 ladder can behave like a 
spin-1 chain, or a spin-1 ladder like a spin-2 chain. 

Assuming that the diagonally situated nearest neighbour spins are also coupled, it was 



shown |13[| that this coupling makes the hidden topological long-range order characteristic 
for the Haldane phase even stronger. 

An alternative way to introduce interchain coupling is to choose the following form for 
the interaction, 



Wi = E n J L( cr t T i+l + °i r , 



T 1 „z^z 



E o J'xyi^i+l + n <i) + J'zTi°l 



+ 1 



(3) 



In this paper we look at the properties of the model defined by Eqs. (|T]) and (|3]). In a special 
case this is equivalent to one of the models studied by White |T5fl , but we will consider a 
larger parameter space, although the anisotropy ratio will be taken to be the same for the 
intrachain and interchain couplings, 

Jz/ Jxy = J z / ' J X yi (4) 

and therefore a single parameter 

J' V 

\ = _pt = l± (5 ) 

" xy <J z 

will be used to characterize the strength of the interchain couplings. 

An interesting feature of our model is that by interchanging the spins between the two 
legs on every even or odd site, the intrachain and interchain couplings change their role. At 
A = 1, where the intrachain and interchain couplings are equal, the model is transformed 
into itself. For other values of A, the energy scale has to be changed. Taking this into 
account, we find that for A > all energy levels satisfy a simple self-duality relationship 



m 



E(X) = XE(1/X). (6) 

This allows to connect the weak- and strong-coupling limits of the model. For A < 
the model would include ferromagnetic and antiferromagnetic interchain couplings, which 
requires the modification of Eq. (||). 



The self-dual A = 1 point plays a special role in other respects as well. At this point, 
with our choice of the couplings, the Hamiltonian of a single chain with composite spin is 
recovered, 

i 

where Si = Oi + fi. It is obvious, that if both a and r are spin-1/2 operators, this model is 
not a true S = 1 Heisenberg model, since the two spins can form not only the triplet, but 



also the singlet combination. It has been shown [14[], however, that the low lying levels of 
the Hamiltonian in Eq. (0) and that of the true 5 = 1 Heisenberg model are identical, and 
therefore they have the same phase diagram. 

The phase diagram expected for the model in the plane spanned by the anisotropy, 
Jzl J-xy and A, is shown schematically in Fig. In drawing the phase boundaries we have 
taken into account that the ferromagnetic phase appears for J z /J xy > 1, independently of 
A. On the other hand the boundary between the Haldane and antiferromagnetic phases is 
expected to depend on A. At A = the ladder model consists of two decoupled spin-1/2 
chains, and therefore the Haldane phase does not exist. At A = 1, according to previous 



finite-size calculations ||i6|-|T9|j on the equivalent 5 = 1 chain, the massive Haldane phase 
exists in a finite range of anisotropy. The numbers obtained for the critical values of the 
anisotropy are somewhat different in the different calculations, but the best estimates give 
the range 

- 1.18 < J z /J xy < 0. (8) 

The value of the critical anisotropy, (J x /J X y)ci — —1-18, where the transition from the 
Haldane phase with singlet ground state to the doubly degenerate antiferromagnetic state 
occurs, is not determined by any symmetry. It is an Ising type transition. Close to it, on 
both sides of the transition point the gap opens linearly. 

It follows from the duality relationship for the energy levels that if A c is a critical point, 
then 1/A C should also be a critical coupling. The phase boundary of the antiferromagnetic 
phase should therefore start from {J z /J xy ) c i = — 1 at A = 0, returning to this value when 
A — > oo, and passing through (J z /J xy ) c i — —1.18 at A = 1. 

Looking at the boundary between the Haldane and planar phases, the massive phase 
appears from the critical planar phase via a Kosterlitz-Thouless type transition. In the 



continuum limit of the 5=1 chain the value (J z / J xy ) C 2 = is obtained |20[ for the critical 
anisotropy. This particular value suggests that it is determined by a hidden symmetry of 
the model, and therefore it is expected that this same critical value could be found in the 



S — 1/2 spin ladder model at any A. This would mean that for — 1 < J z /J X y < the 
interchain coupling is a relevant perturbation. An arbitrarily small coupling will already 
generate a Haldane gap. The corresponding phase boundary is shown in Fig. [I] by a vertical 
straight line at J z /J X y = 0. 

This natural assumption becomes questionable, however, if we compare this phase dia- 
gram with that of the spin-1 ladder. Based on similar considerations we can argue that at 
A = the behavior is that of two decoupled spin-1 chains, while at A = 1 the properties of 
the spin-2 Heisenberg model should be recovered. Since both are integer spin models, the 
Haldane phase exists in both cases in an extended range of anisotropy around the isotropic 
antiferromagnetic point. In the spin-2 model this range is, however, much narrower [^TJ than 



in the spin-1 case, so clearly the boundaries of the Haldane phase to both the planar and 
the antiferromagnetic phases should strongly depend on the interchain coupling, as shown 
in Fig. |. 

Moreover the value of the Haldane gap is much smaller in the S = 2 spin chain than for 
5=1. This is known reliably for the isotropic antiferromagnetic model only [|T^,^], but 
should be true in the anisotropic well. We show schematically in Fig. ^ how the 

antiferromagnetic and Haldane gaps depend on the anisotropy for the spin- 1/2, spin-1 and 
spin-2 Heisenberg models. In all cases the antiferromagnetic state has a doubly degenerate 
ground state with a finite energy gap to the spin wave excitations. For S = 1/2 the transition 
into this state is of Kosterlitz-Thouless type, and the gap opens exponentially slowly. For 
spin-1 and spin-2 chains the transition from the Haldane phase is of second order, so the 
gap opens linearly on both sides. The vanishing of the Haldane gap at (J z /J xy )c2 happens 
again in a Kosterlitz-Thouless-like manner. 

By coupling now two spin-1 chains into a spin-1 ladder, the effect of the interchain 
coupling is expected to be opposite to that discussed above. The Haldane gap decreases or 
it might even vanish when this coupling is switched on. This raises the question whether 
in the spin- 1/2 ladder model the interchain coupling is in fact a relevant perturbation for 
—1 < Jz/Jx y < 0, and perhaps the phase boundary between the Haldane and planar phases 
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does not go in the way discussed above. An alternative possibility is shown in Fig. [I] by a 
dashed line. 

III. NUMERICAL PROCEDURE 

In this paper we try to elucidate this problem by looking at the degeneracy of the ground 
state and the generation or disappearance of gaps when the interchain couplings are switched 
on. In order to do this we have applied the DMRG method [|l^] to the model defined by Eqs. 
(H) and (Q) using the anisotropy J Z J J xy and the relative strength of the interchain coupling, 
A as the two parameters. 

It has been found earlier in various applications, that the DMRG method gives better 
results for systems with open boundary condition (OBC) than for periodic boundary con- 
dition (PBC). On the other hand, when OBC is used, the free spins remaining at the ends 
in the valence-bond like state of the ladder model introduce extra energy-level degeneracies, 
and therefore the analysis of the spectrum becomes somewhat more difficult. In that region 
of the couplings, where the gap is very small, in order to get more accurate results, we have 
used OBC, otherwise we have done calculations with both boundary conditions. 

The major limiting factor in the application of the DMRG method is that only a subset 
of all possible states is kept, and the eigenvalue matrix is diagonalized in this subspace only. 
The accuracy can be improved, if first the Hilbert space is block diagonalized by taking into 
account the natural symmetries of the Hamiltonian, and the DMRG procedure is applied 
to the blocks separately, considering the same fixed number of states in each block. When 
using OBC, translational symmetry is lost, so the total momentum is not a good quantum 
number. In the anisotropic model we also loose the SU(2) symmetry. Therefore only the z 
component of the total Spill, Srp can be used to classify the energy levels. 

It is easily seen that the ladder model Hamiltonian has two further symmetries, two 
mirror planes that are perpendicular to the plane of the ladder. One goes across the middle 
of the rungs of the ladder. It transforms the a and r spins into each other (cr-r symmetry). 
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The other cuts the legs of the ladder in the middle, reflecting the spins on the left onto 
spins on the right (left-right symmetry). Moreover, the = subspace, which will be 
of particular interest for us, has an additional symmetry, the spin reversal or up-down 
symmetry. The corresponding operators will be denoted as V a -T, Vi-r, and V\i- The 
eigenstates of the system can be either symmetric or antisymmetric under these operations. 
An even parity state will be labelled by + while an odd parity state by — . 

It is advantageous to include these symmetries in the DMRG procedure by choosing wave 
functions that are eigenstates of the parity operators. Using the notation of White |23| , the 
superblock wave function tp(aiSi + iSi + 2/3i+3) is formed out of the states ai and $+3 of the 
two blocks and the states si + i and S[ +2 of the two new added spins, 



1>±) 



a l a l + l 



A, 



"1+2^1+3 

x {ip(aisi +1 si +2 (3i+-s) ±Vij(aiSi + iSi +2 Pi +3 )} 



(9) 



.4 



a t Si +1 s l+2 l3i +3 



is a normalization constant which takes the value 1/2 or l/y/2, depending on 



whether the state obtained by the symmetry operation is identical with the original one or 
not. 

The er-r and spin reversal operations are products of local operators that act on the spin 
states of a single rung, so they can be used even if the left and right blocks are not of equal 
length. The four states |||), | ||)> ITI) an d | TT) °f the ° an d t spins on the same rung 
span a basis on which 
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In each iteration cycle of building larger blocks the operators are renormalized as P l+1 = 
OP l O^ and they are stored with the other block operators. The O matrix contains the new 
truncated basis states of the subsequent iteration cycles. A local operator acting on a block 
state gives 
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P l \ai) = a ai \a[) , (11) 

where the coefficient a a[ is 1 or —1 for the symmetric or antisymmetric combination, respec- 
tively 

In the first step of the infinite algorithm, when the block consists of the states of a single 
rung, \cei) — |sj+i) and a ai = 1. The iterated matrix P l+1 has one element in each row and 
that value gives a/ + i. Therefore, the local symmetries simply give 

'Plocali ; ( ( Xl s l+l s l+2pl+3) = 

aaiap l+3 ^( a l s 'l+l s l+2Pl+3)- ( 12 ) 

On the other hand the operator Vj_ r mixes the states of different rungs and can be used 
only when the two blocks are of equal length. When acting on a superblock wave function, 
it gives 

Vi_ r i/}(aiSi + xSi + 2(3i +3 ) = if>{Pi + 3Si +2 si +1 ai). (13) 

The total parity of the wave function is the product of the left-right, f j and cx-r parities. 
In what follows, unless otherwise specified, only the total parity will be given. It is worth 
mentioning that if several symmetries are used, the factor A aiSl+lSl+2 p l+3 has to be included 
for each symmetry. This ensures that every configuration is taken into account only once. 

All these symmetry operations reduce the size of the Hilbert space of the superblock 
Hamiltonian by a factor of about two. This decreases the memory requirements and allows 
to keep more spin configurations, which improves the accuracy of the DMRG method. The 
local symmetries can be included into both the infinite- and finite-system algorithms. The 
left-right symmetry can be used in the infinite-system method only, but by keeping more 
states it would produce better starting vectors |23| for the diagonalizaton procedure in the 
finite-system method. 

Using these symmetry considerations the Hilbert space is split into subspaces character- 
ized by the z component of the total spin and the parity. In order to obtain the ground 
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state and some low lying states of the ladder, in principle one has to determine a few low 
lying eigenstates of all subspaces. For a finite chain with iV sites the n th energy level of the 
Sj. = s sector with parity p will be denoted by E S ^ P (N). Similarly the energy gap between 
the states with energies E s i yn ^ p i and E s ^ p will be denoted by £± s ,n,p;s',n',p>{N). 

Except for the ferromagnetic regime, in the spin-1/2 chain the lowest lying energy level 
belongs to the = sector and it has odd parity under left-right and spin reflection 
symmetries for N mod 4 = while it has even parity for chains with N mod 4 = 2. The 
parity of the first excited state is opposite to that of the ground state and the parity changes 
for every higher lying level. Therefore, in the infinite-lattice algorithm of the DMRG method 
the parity has to be changed in every iteration cycle. 

This difficulty does not arise either in the spin-1 chain or in the ladder models. In 
the spin-1 chain the ground state was found to be in the = sector with even parity 
under left-right reflection for all chain lengths N. In the ladder model beside the left to right 
reflection we also have the cr-r symmetry. The ground state of this model has positive parity 
under both reflections for all non- negative values of A. The excited states also preserves their 
parity as a function of A. Therefore in this case the same symmetry combination has to be 
used in every successive step of the DMRG algorithm. 

Because the energy gaps are expected to be rather small close to the phase boundaries, 
we have used the finite-lattice method version of DMRG with two or three iteration cycles. 
We calculated the energy gaps for ladders with N = 16,32,48,64,82,100,128 sites, and 
used a finite-size scaling procedure to extrapolate to infinite system. Due to our restricted 
computational facility the number of states M representing the block in the DMRG method 
could be chosen between M = 100 and 200 states. The truncation error was worst close 
to the critical points. For the second and higher excited states it was about 10~ 6 — 10~ 7 , 
corresponding to a real error of about 10~ 3 . 

Several formulae have been proposed to extrapolate the finite-size results to the ther- 
modynamic limit. When a finite gap is expected, the N — > oo limit of the gap can be 
obtained by fitting the energy difference of the various E SjUjP energy levels to the form 0] 
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A(iV) = A + A exp(— N/£). On the other hand, for open systems, where the finite size effects 



are more pronounced, the form |22| A(N) = A + A/N or an inverse squared dependence 
|f24| A(N) = A + A/N 2 should be used. The former one is known to give a lower-bound 
estimate of the gap. 

In most cases we have used the formulae with 1/N or 1/N 2 . The exponential dependence 
was assumed when neither the N" 1 nor the N~ 2 fit gave satisfactory result. 



IV. NUMERICAL RESULTS 

In this section we present the results of our numerical calculations. Since the gaps are 
expected to be rather small, resulting in relatively large errors in the extrapolation to infinite 
systems, we are not able to locate exactly the phase boundaries, where the gaps vanish. We 
will therefore restrict ourselves to a qualitative check of how the phase diagram may look 
like. 

Because nothing particularly interesting is expected on the J Z J J xy > side, the system 
is ferromagnetic if J z /J xy > 1 and planar if < J z /J xy < 1 independently of A, in what 
follows we will consider the J z / J xy < region only. We have chosen four different values of 
the anisotropy, namely J z j J xy = —2, —1.1, —0.8, —0.5 and varied A from to 1, to see how 
the character of the ground state changes. The energies will always be measured in units of 

Jxy 



A. The A = and A = 1 lines 

As mentioned above, at A = the properties of two decoupled spin-1/2 chains should be 
recovered, while at A = 1 at least the low energy part of the spectrum is identical to that 
of the spin-1 chain. Our calculations along these lines aimed to test the accuracy of the 
method and to find the relevant subspaces of the Hilbert space. 

Let us consider first the case A = 0. When J z /J xy < — 1, the antiferromagnetic chains 
give a 2 x 2-fold degenerate ground state in the <Sf. = subspace. Accordingly the gap 
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Ao2-,oi+ and two other gaps in the Sj. = sector should go to zero in the limit N — > oo. 
In fact our calculation gives an extrapolated value of the order of 10~ 3 . This shows the real 
accuracy of the numerical calculations. 

The next levels above this fourfold degenerate ground state have one or two spins flipped. 
The corresponding gaps, A n _ )01+ and A 21+>0 i+ were found to be finite, opening slowly as 
the anisotropy increases. Above this gap a continuum of states is found as shown by the 
vanishing of the gap Ai 2+i n_. 

In the planar phase all these gaps should vanish. In fact in the range — 1 < J z /J xy < 
the gaps An_ 01+ and A 2 i +i0 i+ were found to be of the order of 10~ 3 , which we consider to 
be zero. 

Along the line A = 1 the antiferromagnetic phase appears for J z /J xy < —1.18. The 
twofold degeneracy of the ground state is in fact recovered, since the gap A 2-,oi+ scales to 
values of the order of 10 -3 . The antiferromagnetic gap An_ oi+ remains finite. 

In the Haldane phase, —1.18 < J z /J xy < 0, we have again a fourfold degenerate ground 
state, if the ladder has free ends, but here the gaps A 02 _ j01+ and A n _ i01+ should scale to 
zero exponentially. The values obtained for the gap are of the order of 10~ 4 . The next 
levels above them are E 2 \+ and E12+, which cross each other at J z /J xy = —1. Therefore 
the real gap above the ground state is A 2 i+,oi+ for J z /J xy > —1, while it is Ai 2+j0 i+ for 
Jz/Jxy < -1. This gap was found to be 0.147(1), 0.4105(3), 0.223(1), 0.096(3) for J z /J xy = 
— 1.1,-1,-0.8,-0.5. Our result at J z /J xy = —1 agrees quite well with the best estimate 



12] for the Haldane gap of the isotropic S = 1 model. 



B. The antiferromagnetic regime 

Taking now intermediate values for A, we consider first the transition from the Haldane 
phase to the antiferromagnetic phase for J z /J xy < —1. Since the phase boundary is not 
determined by any symmetry, it is expected that for —1.18 < J z /J xy < —1 the antiferro- 
magnetic gap, An_ j0 i+ which is finite at A = 0, should vanish for a finite A c , and then for 
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larger A the Haldane gap Ai 2+j0 i+ should open up. Since the values of the gaps are very 
small, it is difficult to determine accurately the value A c , where this happens. 

As a first attempt, we have investigated the closing of the antiferromagnetic and the 
opening of the Haldane gap as a function of A at J z /J xy = —1.1. The thermodynamic limit 
of the Haldane gap could be determined very well using the N~ 2 dependence for A > 0.1. In 
this region this gap is finite. At the same time the antiferromagnetic gap is found to vanish. 
For A < 0.1, however, no reliable extrapolation procedure was found and therefore we can 
only claim that the transition from the antiferromagnetic to the Haldane phase occurs at a 
very small value of A, between A = and A = 0.05. This indicates, that due to the smallness 
of the antiferromagnetic gap of the spin-1/2 chain, a very small although finite interchain 
coupling is sufficient to drive the system into the Haldane phase. This means that the phase 
boundary starts from the J z /J xy = —1, A = point with an almost horizontal slope. 

Taking a somewhat larger anisotropy, J z /J X y = —1.15, the antiferromagnetic state sur- 
vives to a larger value of A, as shown in Fig. |j. The lower bound estimate of the gap An_ oi+ 
at A = 0.1 is 0.030(3). This point, therefore, belongs to the antiferromagnetic phase. Mov- 
ing up in the phase diagram with A to J z /J xy = —1.15, A = 1 or to the right to the point 
Jz/Jxy = -1-1) A = 0.1, this same antiferromagnetic gap clearly scales to zero, indicating 
that these points are already in the Haldane phase. The phase boundary should therefore 
have the shape indicated in Fig. |l|. 

C. The planar phase 

Next we considered the stability of the planar phase in the — 1 < J z j J xy < regime. We 
have measured the gap A2i+,oi+, which should vanish at A = 0, as a function of A for several 
values of the anisotropy. 

The accuracy of our calculations was checked by comparing the energy values calculated 
for A < 1 and A > 1, using the self-duality relationship of Eq. (|6]). We have got the same 
energies to at least 6 or 8 digits. We show in Fig. [5] the results obtained for J z /J xy = —0.8. 
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At A = the gap scales to 0.008(2) as 1/N. At A = 0.2 an excellent fit can be obtained by 
the 1/N 2 fit for N = 128, 100, 82. The extrapolated gap was found to be 0.057(2). 

The finite- size scaling procedure is not that clear for A = 0.1. Neither a 1/N nor a 1/N 2 
fit is satisfactory. This indicates that the available chain lengths are still not long enough to 
determine A accurately. The estimated value of the gap is between 0.025(3) and 0.035(3). 
Comparing this with the gap at A = 0.2, we could conclude that in the thermodynamic limit 
the gap opens linearly with A c = 0, in the same way |1| as for J z / J xy = — 1. 

For J z / J xy closer to zero, the Haldane gap is too small even at A = 1. At J z / J xy = —0.5 
the extrapolated value of the gap is about 0.096(3). For small values of A of the order of 
0.1 the error of the DMRG method becomes comparable to the value of the gap and no 
reliable extrapolation procedure was found. The results are, however, in agreement with the 
assumption that A c = 0. 

These results indicate that the interchain coupling is relevant in the whole —1 < J z / J xy < 
region. An arbitrarily weak interchain coupling of this type would drive the spin-1/2 ladder 
into the Haldane phase. 



D. The spin-1 ladder model 

The spin-1 ladder model is constructed in the same way as the spin-1/2 model except 
that the Oi and r, operators in Eq. (H) and @ are now spin-1 operators. In the A = limit, 
which in this case corresponds to two decoupled spin-1 chains, the Haldane phase exists for 
a wide range of anisotropy, namely for —1.18 < J z /J xy < 0. At A = 1 this model should 
behave like a spin-2 chain, in which the Haldane phase is squeezed into a very narrow region 



around the isotropic point. The gap itself is also small. According to the best estimate [22] 
it is A = 0.085(5) at the isotropic point. 

Therefore, as shown if Fig. 0, the Haldane gap is expected to vanish at a finite value of 
the interchain coupling when J z /J xy is not too close to —1. In order to see this, we have 
looked at the gap A 2 i +> oi+, which is finite in the Haldane phase, but vanishes in the planar 
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one. It was found that at J z /J X y = —0.7, A = 0.5 this gap was already zero. Although 
we could not locate the boundary precisely, it is clear that, taking into account again the 
self-duality relationship, the phase boundary should be perpendicular to the anisotropy axis 
at A = 1 and have the shape shown in Fig. |^. 

V. DISCUSSION 

In this paper we have studied the stability of the Haldane phase in a magnetic ladder 
model. Our special choice of the coupling between the chains allowed us to investigate and 
smoothly interpolate between both the half-odd integer and integer spin Heisenberg models. 
We have applied the DMRG method which was improved to reduce the computational time 
and to increase the accuracy. 

We have found that for the 5 = 1/2 + 1/2 composite spin model the phase boundary 
between the antiferromagnetic and the Haldane phase is a smooth curve and the value of A c 
is a function of the anisotropy parameter, J z /J xy . For the boundary between the Haldane 
and planar phases our results indicate that it is at J z /J xy = independently of A. For 
< Jz/Jxy < 1 the system remains in the planar phase, while for —1 < J z /J xy < the 
interchain coupling is relevant. The critical value where the Haldane gap starts to open up 
is A c = 0. 

This is, however, not the only possibility. Our results are not in contradiction with the 
assumption that A c is finite, although small, except when J z /J xy is rather close to 0, as 
shown by the dashed line in Fig. [l| 

We should also mention that since the point A = 1, J z /J xy = on the boundary between 
the Haldane and planar phases was obtained in the continuum limit, it cannot be excluded 
that in the lattice model this boundary at A = 1 is not exactly at J z j J xy = 0, but at a small 
negative value. In this case the boundary between the Haldane and planar phases should 
be slightly deformed from that shown in Fig. [I]. 

On the basis of our calculations, unfortunately, we are unable to decide, which scenario 
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takes place. The argument that gave (J z /Jxy)c2 = at A = 1 would give a phase boundary 
which is independent of A, meaning that the interchain coupling is relevant for any J z j J xy 
in the range — 1 < J z /J xy < 0. The same reasoning might then indicate that the interchain 
coupling is relevant also in the case when four S — 1/2 chains are coupled. With an 
appropriate choice of the couplings this would, however, generate a four-leg ladder that 
behaves like an S = 2 chain. In this model the phase boundary between the planar and 
Haldane phases cannot be independent of A, in the same way as in the two-leg S = 1 ladder. 
Further calculations, keeping probably close to 1000 states in the DMRG procedure, are 
necessary to decide between these possible scenarios and to determine the precise shape of 
the phase boundaries. 
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FIGURES 

FIG. 1. Schematic phase diagram of the S = 1/2 ladder model in the (X,J z /J xy ) plane. The 
dashed line shows an alternative phase boundary between the Haldane and planar phases, as 
discussed in the text. 

FIG. 2. Schematic phase diagram of the S = 1 ladder model in the (A, J z j ' J X y) plane. 
FIG. 3. The energy gap of the S = 1/2, 1 and 2 Heisenberg chains as a function of the anisotropy 

Jz/Jxy- 

FIG. 4. The antiferromagnetic gap An_ ) oi+ as a function of iV -1 at (a) J z /J xy = —1.15, 
A = 0.1 and A = 1 and (b) J z /J xy = -1.1, A = 0.05 and A = 0.1. 

FIG. 5. The energy gap A2i+,oi+ as a function of N~ 2 at A = 0.2, J z /J xy = —0.8. 
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